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ABSTRACT 



Aims. We revisit particle acceleration in the polar cap region of a neutron star by taking into account both general relativistic effects 
and the presence of toroidal oscillations at the star surface. In particular, we address the question of whether toroidal oscillations at 
the stellar surface can affect the acceleration properties in the polar cap. 

Methods. We solve numerically the relativistic electrodynamics equations in the stationary regime, focusing on the computation of 
the Lorentz factor of a space-charge-limited electron flow accelerated in the polar cap region of a rotating and oscillating pulsar. 
To this extent, we adopt the correct expression of the general relativistic Goldreich- Julian charge density in the presence of toroidal 
oscillations. 

Results. Depending on the ratio of the actual charge density of the pulsar magnetosphere to the Goldreich- Julian charge density, we 
distinguish two different regimes of the Lorentz factor of the particle flow, namely an oscillatory regime produced for sub-GJ current 
density configurations, which does not produce an efficient acceleration, and a true accelerating regime for super-GJ current density 
configurations. We find that star oscillations may be responsible for a significant asymmetry in the pulse profile that depends on the 
orientation of the oscillations with respect to the pulsar magnetic field. In particular, significant enhancements of the Lorentz factor 
are produced by stellar oscillations in the super-GJ current density regime. 

Key words. Magnetic fields, Relativistic processes. Stars: neutron, (Stars): pulsars: general 



1. INTRODUCTION 

Understanding the physical mechanism behind particle 
acceleration, and hence electromagnetic emission, in the 
magnetospheres of neutron stars remains an unsolved topic 
of pulsar physics, more than forty years after their discov- 
ery. A completely force-free (FF) magnetosphere with no 
parallel component of the electric field along the magnetic 
field lines, clearly does not allow for the acceleration of par- 
ticles, and local viola tions to the FF condition are therefore 
exp ected ( Kalapothar akos et al.l 1201 iT) . This was first noticed 
by ISturrocki (il971i) . who proposed that the acceleration of 
particles occurs in a small vacuum gap, where E ■ B + 0, 
that is confined to the polar region above the star surface, the 
so-called polar cap, whose edge is defined by the locus of the 
last closed magnetic field line. This seminal idea was developed 
and investigated by several author s ove r the y ears, both in 
flat s pacetime jRuderman & Suther land 19751; Fawlev et al 



19771: 'Scharlemann et al. 1978; "Xrons & Scharlemann' 



1979 : iCheng & Rudermaa ,1980; .Daughertv &l Hardingj 
19821 IShibatal Il997r and i n a gen er al relativistic frame 



work (LMuslimov & Tsygan" 1992; "iVIuslimov & Hardin 



19971; iHarding & Muslimov 
Morozova et al.ll2008l) . 



1998; Sakai & Shibata ,20031 



A relevant result obtained by these works is that the parti- 
cle motion can show an oscillatory behavior depending on the 
charge current density. This behavior has been indirectly con- 
firmed by numerical simulations. By performing time-dependent 



particle -in-cell simulati ons of electrons extracted from the star 
surface, lBarzilavl(l20Tl found that the Lorentz factor of the ac- 
celerated electrons reaches the analytical estimates of 10^ and 
follows an oscillatory pattern that increases with the distance 
from the stellar surface. Although these numerical analyses are 
very promising, numerical simulations of neutron star magne- 
tospheres have not yet entered a mature enough stage where 
physical mechanisms such as plasma effects, non-ideal magne- 
tohydrodynamics, pair-creation, and radiation processes can be 
properly taken into account. As a result, analytical investigations 
continue to play a crucial role in elucidating the physics of pulsar 
magnetospheres. 

Among the physical effects deserving particular attention is 
the possibility that neutron star oscillations, most likely excited 
during a glitch phenomenon (i.e. a sudden change in the rota- 
tional period), propagate into the magnetosphere, thus affecting 
the acceleration properties in the polar cap region. The idea that 
neutron star oscillations may induce high energy emission in 
neutron star magnetospheres w as developed in a series of papers 
by Timokhin and collaborators dTimokhin et al.ll2000t 'Timokhinl 
2007; Timokhin et al. 2008), after the first pioneering in vestiga- 
tions by McDermott et al. ( 198?), iMuslimov & TsvganI (Il986l) . 
and iRezzolla & Ahmedov. (.20041) . who considered the case of 
an oscillating neutron star in a vacuum. A few years ago we ini- 
tiated a program aimed at assessing the impact of toroidal os- 
cillations of the star surface on the properties of the external 
magnetosphere, with genera l relativistic effects proper ly taken 
into account. Starting from lAbdikamalov et aP (l2009l) . where 
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th e theoretical basis of ou r approach was developed, we showed 
in lMorozova et al that the electromagnetic energy losses 

from the polar cap region of a rotating neutron star can be sig- 
nificantly enhanced if oscillations are also present, and, for the 
mode {l,m) = (2, 1), these electromagnetic losses turned out to 
be a factor ~ 8 larger than the rotational energy losses, even for 
a velocity oscillation amplitude at the star surface as small as 
7/ = 0.05 D. R, where Q and R are the an gular velocity and the ra - 
dius of the neutron star, respectively. In lMorozova et al.l (1201 2h . 
on the other hand, we considered the conditions for radio emis- 
sion in magnetars and we found that, when oscillations of the 
magnetar are taken into account, radio emission from the mag- 
netosphere is generally favored. The major effect of oscillations 
is to amplify the scalar potential in the polar cap region of the 
magnetar magnetosphere, which implies that the death-line in 
the P - P diagram shifts downwarqj 

As a natural follow-up of our investigations, we address in 
this paper the question of whether toroidal oscillations at the star 
surface can affect the acceleration prope rties in the polar cap re- 
gion. We rely on our results presented in lMorozova et al.l (20101) 
revisiting some of the conclusions found bv ISakai & Shibatal 
(l2003h . In particular, we expect that an oscillating magneto- 
sphere can alter the parallel component of the electric field that is 
responsible for particle acceleration just above the star surface. 
The size of oscillation amplitude necessary to produce signifi- 
cant p hysical effects is still a matter of debate. Ti mokhin et al.. 
(I2008h . for example, argue that noticeable effects in the spectra 
of soft gamma repeater giant flares can be produced when the 
amplitude of the neutron star surface oscillations is of the order 
of ~ 1 % of the star radius. As we show in the paper, on the other 
hand, significant effects on the Lorentz factor are visible even 
with (normalized) velocity oscillation amplitudes K - fj/QM as 
low as 0.02. Since the range of oscillation freque ncies is quite 
broad, typically between 18 Hz and 1800 Hz (Ste iner & WattsI 
l2009l) . such low values of K can be obtained, for example, by 
adopting an oscillation frequency that is 20 times higher than 
the stellar rotation frequency Q., and a spatial oscillation ampli- 
tude that is only 0.0017?. We anticipate, therefore, that the effects 
of an oscillating magnetosphere on the acceleration of particles 
close to the stellar surface become important, even for very small 
oscillation amplitudes. 

The plan of the paper is the following. In Sec.|2] we describe 
the basic features of the physical model, while Sec.|3]is devoted 
to the presentation and discussion of the results. Sec. |4] contains 
the conclusions of our work. In addition to the Gauss units for 
the magnetic field, we set c = G = 1. In this way, any quan- 
tity, including the magnetic field, is geometrized, namely there 
is only one unit of measure, which is the gravitational radius 
rg - GM/c'. We keep c and G in an explicit form in those ex- 
pressions of particular physical interest. Appendix |A] describes 
the extended geometrized system of units adopted in the numer- 
ical solution of Eqs. (|8]l and (|9]l. 



2. THE PHYSICAL MODEL 

2. 1 . Particle dynamics in tine pulsar polar cap 

We assume that the space-time outside a representative neu- 
tron star of mass M, radius R, and angular velocity is 



' However, we also found that when the compactness of the neutron 
star is increased, the death line shifts upwards in the P - P diagram, 
pushing the magnetar into the radio-quiet region and explaining the ob- 
servational evidence that most of the known magnetars are radio-quiet. 



given by the Hartle & Th orne metric in the slow-rotation 
hmit (iHartle &Thornell968h . i.e. 

ds^ = -N'^di^ + N'^dr^ - Icou^r^ sin^ 6dcf>dt + 



r^ide^ + sin^ 0d(p~) , 



(1) 



where - yjl - is the lapse function, wlt = is 

the Lense-ThiiTing angular velocity, and J - /Q is the to- 
tal angular momentum of the neutron star measured from in- 
finity, with I being the moment of inertia. The frozen dipole- 
like magnetic field of the pulsar has the following orthonor- 
mal co mponents in spherical coordinates (r, 0, 0) (see, for in- 
stance, Ginzbiirg & Ozernoilll965l: IWasserman & Shapiroll 19831 : 



iRezzolla et al.ll200ll) 



Br 



fio^rr- — COS 6, 



:BoN 



where 



7(1) 



Inl 



(l-f)/(l) 



f + f 1 + A 



■ sin 9, 



(2) 



(3) 



(4) 



In the expressions above, 77 = r/R is the normalized radius, s = 
2M/R is the compactness parameter of the neutron star, and Bo is 
the value of the magnetic field at the pole of the star, which is of 
the order of lO'^G for a typical pulsar Ove r the years, a series 
of fu n damental works (see, am o ng oth ers, Goldreich & Julian 
1969 : Ruderman & Sutherland! 1975 : Arons & Scharlemanii 



T979[ iMestel et a l.l 119851: iFitzpalrick & Mestell [19881: iMichei 
Muslimov &Harding|ll997l: iBeskin 2oM^ 



1 99 It iMuslimov & Harding|II997l: [Beskin 2010|) have helped 
to reveal the mechanism behind the pulsar magnetosphere 
generation. A strong surface magnetic field and a relatively 
rapidly rotating pulsar lead to the generation of an electric field 
in the vicinity of the pulsar. In the framework of a space- 
charge-limited-flow accelerator that we assume here, the par- 
allel co mponent of t he electric field £[1 at the stellar surface 
is zero (lBeckeijl2009l) . Therefore, if the surface temperature is 
higher than the thermionic emission temperature, particles can 
leave the star, but with non-relativistic speeds. Once above the 
star surface, electrons are accelerated to ultra-relativistic veloci- 
ties on a short distance scale by the non-zero Zsy. This first gen- 
eration of electrons, called primary electrons, emit curvature y- 
quanta, which, in turn, create electron-positron pairs, initiating 
the cascade process that forms the pulsar magnetosphere. 

The charge density of the pulsar magnetosphere, which is 
required in order to completely screen the electric field com- 
ponent parallel to the magnetic field of the pulsar, is called the 
Goldreich- Julian charge density. The actual charge density p of 
the pulsar magnetosphere slightly deviates from the Goldreich- 
Julian charge density and the deviati on can be rou ghly estimated 
from the one-dimensional equation (lBeskinll20IOh 

dh 



- P- PGJ,rot , 



(5) 



where h is the distance from the star surface. The expression for 
the Goldreich- Julian charge density of the pulsar, corrected for 
general relativistic effects, was obtained bv lMuslimov & TsvganI 
(II992h and is given by 



PGJ.rot - 



1 m 



Inc Nif /(I) 



co&X + 2^^^^^ sin;^f cos (p 



,(6) 
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where 



Hirf) 



77- 



+ 1 



3e K 
2t] 2rf 



1 



(7) 



while X is the inclination angle between the axis of rotation of 
the pulsar and the magnetic moment of the dipole. Moreover, 
K = e/3 - (a>ur)*/^, where the asterisk denotes a value at the 
surface of the star, and /3 - I/Io is the normalized moment of 
inertia, where /o = MR~. Since the polar cap is located in the 
vicinity of the magnetic pole, in Eq. ^ the polar axis of the 
spherical coordinates (r, 6, (p) is chosen to be directed along the 
magnetic dipole axis, which also explains why the angle x does 
not appear in the expressions for the magnetic field components 
given in Eqs. (|2]i and ([3]). 

Within this framework, it is possible to solve the equations of 
motion of the electron accelerated in the vicinity of the polar cap. 
Under the ass umption t hat th e cha rge flow reaches a statio nary 
configuration, IShibatal ( Il997h and ISakai & Shibatal ( l2003h for- 
mulated this problem in the form 



d4>\ r(i' + 1) 



ds 



Nd_ 



— (N \ - 



Ns^ 



B 



where 



(8) 
(9) 

(10) 



is a normalized radial coordinate, while v is the velocity of the 
electron. In equations ([8]) and (|9]l, additional normalizations are 
performed as 



2kN(s)J(s) 



4.is) = -0(i) 

m 



2nN (s) c pG]jot(s) ^ 

QB(i) 

' k\ 3 
1 -\ cosx H — H{rf)G s\nx cos < 



(11) 
(12) 

(13) 



where J is the current density in the pulsar magnetosphere, <b is 
the scalar electric potential, m and e are the mass and the charge 
of the electron, respectively, y is the Lorentz factor of the ac- 
celerated electrons and B s |B| = (B? + Biyi- ^ B;- + 0(6^). 

^ 

Moreover, we stress that, in order to derive the equation (|8]l from 
the relativistic Poisson equation, the scalar potential <1) has been 
expanded in the angular directions as 



1> = 2<1>W(0,0), 



(14) 



l',m' 



where the spherical orthonormal functions Yi„,{6, <p) are the 
eigenfunctions of the Laplacian operator in spherical coordi- 
nates. Only modes of the polar cap scale have been considered, 
with r » 7r/0o, where 



©0 = arcsin 



R 



RLcf(l) 



1/2 



(15) 



is the polar ang le of the last open magnetic field line at the sur- 
face of the star dMuslimov & TsvganI 19921) . while Ric - c/Q.is 



the radius of the light cylinder In our computations, for typical 
values of ©o in the range ^ [2°, 3°], we have /' * [60, 90]. The 
polar trajectory of the accelerated particles is described by the 
curve 



(16) 



and remains b ounded by the polar an gle of the last open mag- 
netic field line dMorozova et al.ll2010h . i.e. 



0(s) < = arcsin < 



/(I) 



sin ©n 



(17) 



2.2. The effects of pulsar oscillations 

In our analysis, we also consider the case when the rotating pul- 
sar is subject to toro idal oscillations, with the orthonormal ve- 
locity field given by dUnno et al.ll 19891) 



6v' = .^0 



sin 6 



d^YUO,'l>),-dgYUe,'f>)\fi(r)e 



(18) 



where oj is the real part of the oscillation frequency and fj is the 
radial eigenfunction expressing the amplitude of the oscillations. 
We assume that the oscillations of the stellar crust can induce 
oscillations in the plasma magnetosphere, at least in the close 
vicinity of the surface, which is our main region of interest. As 
a result, in the rest of our analysis we neglect the radial depen- 
dence of the oscillation amplitude, i.e. we assume fjir) = rjiR). 

In principle, an additional degree of freedom should be in- 
troduced, related to the possible uncertainty in the orientation of 
the oscillations with respect to both the axis of rotation and the 
pulsar magnetic field. However, we adopted here a simplified ap- 
proach and, as implied by Eq. ( fTSl l, oriented the oscillation mode 
axis along the z- axis of the magnetic dipole, while the rotational 
axis is inclined. 

The Goldreich-Julian charge density of a ro t ating and os- 
cillating_pulsar was obtained by iMorozova et al.l (1201 Oh and is 
given by 



PGJ - PGJ.rot + PGJ.osc - 

1 1 Boe-'"' 



PGJ.rot 



4-ncRT]'^ @\rt) N f(l) 



(19) 



After using the expression above, one can easily obtain the cur- 
rent j in the presence of both rotation and oscillations as 

= |l - ^jcos;r + ^//(77)0sin^cos0 + 

1 / f(n^\^ ff" 

- — —^Kri'^^-Hil + \)Ai,„ , (20) 

2 1/(1)/ @l ' ' ' 

where the amplitude of the oscillation has been parametrized in 
terms of the small number K = fi{R)/QR, giving the ratio of the 
velocity of oscillations to the linear rotational velocity of neutron 
star. For any particular mode (/, m), we used the approximation 
Yii„{6, <p) Si Aii„{<p)ff^\ which is valid in the limit of small polar 



^ We were careful not to confuse /' and in' (which are used in the 
expansion of the scalar potential in Eq. jl4b). with / and in used in the 
expansion of the perturbation in Eq. dlSb . 
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angles G, where the terms Ai,„{(f)) appearing in Eq. ( l20t have real 
parts given by 



(21) 
(22) 
(23) 
(24) 
(25) 




3./ — cos< 



We note that in all our analysis we only consider constant 
ratios j/ y,, namely we assume that the ratio of actual charge 
density to the Goldreich- Julian charge density is constant. This 
should be regarded as a very good approximation as long as the 
acceleration mechanism is confined to a short distance from the 
stellar surface. 

We stress that the two most plaus ible mo d els of pul- 
sar magnetosphere, n amely those of iMestell (Il981h and 
lArons & ScharlemannI (Il979l) . predict different distributions of 
the charge density in the polar cap region, and hence different 
conditions for pair creation. In the model of Mestel, a super-GJ 
charge density is formed along the magnetic field lines that curve 
away from the rotational axis, while along the field lines curving 
towards the rotational axis a perfect screening of the accelerating 
electric field component can be achieved and the charge density 
is equal to the Goldreich-Julian charge density. In contrast, in the 
model of Arons, the region where magnetic field lines curve to- 
wards the rotational axis has a super-GJ charge density near the 
star surface, and a sub-GJ charge density away from iQ One of 
the main differences between the two models, which is responsi- 
ble for different predictions about the charge density in the pulsar 
magnetosphere, lies in the boundary conditions of the problem: 
in the model of Mestel, in particular, the charge density is taken 
to be equal to the Goldreich-Julian charge density at the surface 
of the star, while Arons suggests that E\i = far away from the 
surface. 

However, the present lack of a fully self-consistent model of 
pulsar magnetospheres motivates a pragmatic approach, accord- 
ing to which the charge density can be treated as a free parame- 
ter. 



3. RESULTS 

In a first series of calculations, we have solved the system of 
equations given by Eqs. ([S])-® for the case of a rotating but 
non-oscillating pulsar. The two equations Eqs. ([8]) and ^ have 
been reduced to a system of three ODEs in the unknowns (p, 
dip Ids, and Ny, which has been solved with a standard fourth- 
order Runge-Kutta method. The main parameters of our fiducial 
pulsar model are 7? = 10 km, M = I.4M0, Q. = 2;r/(0.1s), Bo = 
1.0 X lO'-G, K = 0.1, and x = 30°. We note that, by selecting 
K = 0.1, the moment of inertia of the pulsar is given by ///q ~ 
0.24. The initial value of the Lorentz factor at the star surface 
is 7, = 1.01, while the initial 0* is bounded by ©o, namely is 
0<e,< 00- 



108 ^ 



105 ^ 



10" r / 



j-0,98j, n 
J = 99L 
j=1.0L 




0.002 0.004 

r/K- 1 



0.006 



Fig. 1. Dependence of the Lorentz factor on the ratio jl y, for a 
neutron star with M = I.4M0, R ^ 10 km, P = Q.\s,x = 30°, 
Bo= l.Ox 10'"G, = 0°, 00 = 2°, and r, = 1.01. 



70 

60 H 
50 ^ J - h 



B„= l.Ox 10" G ^ 
B„=l-OxlOi" G I 
B„=1.0xlO" GH 




0.001 



o.ooa 0.003 

r/R- 1 



0.004 



' On the other hand, this type of solution does not apply in the region 
where magnetic field lines curve away from the rotational axis. 



Fig. 2. Lorentz factor dependence on the intensity of the mag- 
netic field for a neutron star with M - 1.4Mo, R - \Q km, 
p = 0.\s,x = 30°, 6, = 0°, ©0 = 2°, = 1.01. Top panel: 
i = 0.98 Bottom panel: ; = 1.01 7.. 



Fig.[T]reports the Lorentz factor as a function of the normal- 
ized radial coordinate rjR- 1 for differe nt values of the ratio /'/ /„ . 
Fig.[T]essentially confirms the resul ts of Sakai & Shibatal (l 2003h . 
which were originally predicted bv lMestel et al. " (n'985h . namelv 
the existence of two different acceleration regimes. When j < j, 
small values of the velocity v can make the right-hand side of 
Eq. ([8]) positive, thus producing an increasing potential and an 
increasing Lorentz factor As soon as the velocity reaches a crit- 
ical value, given by the initial ratio jj j, the right-hand side of 
Eq. ^ becomes negative and the potential, as well as the veloc- 
ity, starts to decrease. As a result, the case j < j is responsible for 
the occurrence of an oscillatory regime, with the oscillation fre- 
quency (in the spatial domain) that is strongly dependent on the 
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Fig. 4. Lorentz factor dependence on the normalized amplitude of the stellar oscillations K for the mode of oscillations (/, m) = (1, 1) 
with 0, = 2°, ©0 = 3°,y, = 1.015, and Bq - 1.0 x lO'^G. The two top panels correspond to the case j - 0.98 J„ while the two 
bottom panels correspond to the case j - 1 .001 y,. The left panels show the solution for = 0, the right panels for (p - n. 



ratio il 7/. As shown bv lMestel et"al1 ( Il985h and lShibatal ( Il997h . 
this situation corresponds to the case where the effective charge 
density p-pcj.iot changes sign. Whenp-pGj,iot > 0, the electron 
flow decelerates to increase the charge density; in contrast, when 
P - PGj.rot < 0, the flow accelerates to reduce the charge density. 
A rather different behavior is instead encountered when j > j. 
In this case, the effective charge density is always negative, the 
right hand side of Eq. ^ is always positive, and the potential as 
well as the velocity can only grow. 



ISakai & Shibatal ( l2003h erroneously attributed the abrupt in- 
terruption of the oscillatory pattern, which is manifested by some 



models (see the bottom panel of their Fig. 1), to a purely gen- 
eral relativistic effect. In contrast and as we have verified, this 
effect is purely numerical and caused by v~ possibly becom- 
ing negative during the integration of the system of equations 
Eqs. (O-©. For instance, this is also true for the long-dashed 
blue curve in Fig.[T] where the Lorentz factor decays to unity at 
r/R - 1 » 0.003 after having reached a value of ~ 4100. 
The pathology of the ODE system that we have demonstrated 
could not be solved even when resorting to stiff solvers. As a re- 
sult, the behavior of the Lorentz factor after v~ becomes negative 
cannot be trusted and is not reported in our figures. We also note 
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Fig. 5. Lorentz factor dependence on the normalized amplitude of the stellar oscillations K for the mode of oscillations (l, m) - (2, 1) 
with 0, = 2°, 00 - 3°,7, - 1.015, and Bq - l.Ox lO'^G. The two panels correspond to the case j = 1.001 y',. The left panel shows 
the solution for = 0, and the right panel for cp - n. 




o.ooa 



0.004 



0.006 



r/R 1 



Fig. 3. Lorentz factor dependence on the incUnation angle x for 
a neutron star with M - IAMq, R - 10 km, and P - OAs, 
j = l.OlJ,, e, = 0°, 00 = 2°, J, = 1.01, and Bo = 1-0 x lO'^G. 
The Lorentz factor decreases for larger inclination angles. 



that a value of r/R -I - 0.006 corresponds to a distance of 60 m 
above the star surface, which is smaller than the radius of the 
p olar cap ^ 0o ~ 350 m as in ferred by the polar cap model 
of lRuderman & Sutherla"ndl(ll975h . 

On the other hand, Fig|2] illustrates the dependence of the 
Lorentz factor of the accelerated particles on the intensity of 
the magnetic field. We note that both the accelerating electric 
field £|| and the difference p - pgj.i oi are directly proportional t o 
the magnetic field (See Eq. (56) bv lMuslimov & Tsvganl[T992l) . 
Therefore, if j > jt, corresponding to the case when the eff'ective 
charge density is always negative, the Lorentz factor increases 
with the magnetic field, as shown in the bottom panel of Fig|2l 
On the other hand, if j < j,, corresponding to the case when the 
effective charge density changes sign, the charge density jump 
Ip -pGj.rotl increases with the magnetic field, producing a pattern 
with a higher frequency of oscillation (See Eq. (23) bv lShibatal 
Il997h . This effect is reported in the top panel of Fig|2l 



The dependence of the Lorentz factor on the inclination an- 
gle is shown in Figl3] When j < jt, namely when the oscil- 
latory pattern is produced, the Lorentz factor is insignificantly 
affected by variations inx and the results are not reported here. 
In contrast, when j > jt, the results are more interesting and 
show that smaller Lorentz factors are obtained for larger inclina- 
tion angles although this effect becomes significant only for 
very inclined configurations. 

In a second series of calculations, we have considered the 
effects of stellar oscillations on the Lorentz-factor of the accel- 
erated particles by adopting the prescription given in Eq. ( l20t . 
Figures |4] and |5] show the solutions for the modes of oscilla- 
tions {l,m) = (1, 1) and (2, 1), respectively. As is apparent from 
the two top panels of Fig. H] the presence of stellar oscillations 
does not affect the amplitude of the Lorentz factor of the par- 
ticle for j - 0.98 J. and (l,m) - (1, 1). On the other hand, the 
frequency of the velocity oscillation is smaller for larger K at 
<p -0° (left-top panel), while it is larger for larger K alcp - 180° 
(right-top-panel). As long as the velocity of oscillations for the 
considered modes is proportional to cos (p, the effect will be the 
opposite when cos (p changes sign. The same qualitative behav- 
ior was also encountered for the mode (l,m) = (2, 1), and we 
conclude that it is typical of the case j/ y* < 1 . 

The case j/ jt > 1 is reported in the two bottom panels of 
Fig.|4]for the mode {l,m) = (1, 1) and in Fig.[5]for the mode 
(l,m) = (2, 1), and illustrate that star oscillations can signifi- 
cantly affect the Lorentz factor. We note that the influence of the 
stellar oscillations on the acceleration of the particle strongly 
depends on the azimuthal angle (p. In particular, owing to the 
modulation produced by the term cos (p, larger oscillation ampli- 
tudes K become responsible for a net decelerating effect at an- 
gles (p ~ 0, while they produce a net accelerating effect at angles 
(p ~ 7T. Moreover, a strong dependence on the ratio j/ y, is also 
encountered, as is clearly evident from the two panels of Fig.|6] 
which show the modulation of the Lorentz factor with the an- 
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gle (p for the mode {I, m) - (1, 1) in a representative model with 
K — 0.02. The left panel refers to the case jl /', - 1.001, while 
the rig ht panel has jlj^ - 1 .01 . Interestingly. iMorozova et alJ 
(1201 oh found that the oscillation modes with m — \ considerably 
increase the electromagnetic energy losses from the polar cap 
region of the neutron star, losses that can be several times larger 
than in the case when no oscillations are present. A more effi- 
cient particle acceleration typically corresponds to higher energy 
losses from the neutron star magnetosphere ( Osmanov & Rieger 
I2009h . As a result, we conclude that, at least for the modes 
{l,m) = (1, 1) and {l,m) - (2, 1) that we have considered here, 
larger oscillation amplitudes at the stellar surface may be re- 
sponsible for both stronger accelerations and higher luminosi- 
ties, provided that the current density in the magnetosphere is 
higher than the Goldreich-Julian current density. 

At least in the case of small inclination angles, a noticeable 
influence of stellar oscillations on the conditions for particle ac- 
celeration can be intuitively understood as follows. The acceler- 
ation of charged particles extracted from the surface of the pulsar 
in the polar cap region is determined by the electric field, which 
is generated by the combined effect of strong magnetic fields 
and the motion of the star surface. When stellar oscillations are 
added to the pure rigid rotation, a given velocity distribution at 
the star surface is produced. The rotation velocity in the polar 
cap region is proportional to the small polar angle 0, is exactly 
equal to zero on the axis of rotation, and increases with dis- 
tance away from the axis. In contrast, the velocity of oscillations, 
which is given by Eq. ( fTSl ). in the case m - 1 does not depend 
on the small polar angle and is effectively constant across the 
polar cap region of the star This explains why for these modes 
even small oscillation velocities may have a noticeable effect on 
the accelerating component of the electric field. 



4. CONCLUSIONS 

By numerically solving the relativistic electrodynamics equa- 
tions in the stationary regime of a space-charge-limited-flow ac- 
celerator, we have computed the Lorentz factor of electrons ac- 
celerated in the polar cap region of a rotating and oscillating 
pulsar magnetosp here. Our results confirm some of the funda- 
mental findings of lSakai &~Sh ibata (2003), namely the existence 
of two different regimes: an oscillatory one, which is produced 
for sub-GJ current density configurations, (j/jt < 1) and does 



not produce an efficient acceleration, and a truly accelerating 
one, which is produced for super-GJ current density configura- 
tions (j/ > 1). We have also clarified the numerical origin of 
a stopping effect on the velocity of the electron, which has noth- 
ing to do with general relativistic effects. Stellar oscillations, on 
the other hand, can affect both the absolute value of the Lorentz 
factor gained by the electron and the frequency of the oscilla- 
tory patterns, when the latter is present. Owing to the modula- 
tion produced by the term cos in the equations, where is the 
azimuthal angle, the way in which each mode of oscillation af- 
fects the dynamics is different at different azimuthal positions on 
the polar cap surface, and, in particular, it follows the periodic- 
ity of cos (/). These results are overall consistent with those ob- 
tained by Morozovaet al. (2010), who showed that oscillations 
are responsible for an extra term in the total energy losses from 
the system, because the electromagnetic energy losses are deter- 
mined by the integrated absolute value of the current, flowing 
in the magnetosphere. At least for the modes (/, m) = (1,1) and 
{l,m) - (2, 1), larger oscillation amplitudes at the star surface 
may produce both higher accelerations and higher luminosities, 
provided that the magnetosphere has a super-GJ current density 
configuration. However, these results also suggest that pulsar os- 
cillations may become responsible for a significant asymmetry 
in the pulse profile, which will depend on the orientation of the 
oscillations with respect to the pulsar magnetic field. 

Finally, we stress that some recent investigations have tried 
to connect the models of stellar oscilla tions with the obser- 
vational data available f or pulsars (iRosen & Clemeni 120081 : 
iRosen & Demoresil201 lb iRosen et al.ll201 ih . The scenario that 
is emerging from these studies is that the presence of stellar os- 
cillations creates different kinds of "noise" in the clock-like pic- 
ture of pulses, i.e. changes in the pulse shape, changes in the 
spin-down rate, and the switching between different regimes of 
pulsar emission. Ivan Leeuwen & TimokhinI (l2012h . in particular, 
revisited the issue of drifting subpulses in the pulsar profiles, 
reaching the important conclusion that the rate of subpulse drift 
is proportional to the latitudinal derivative of the scalar potential 
across the pulsar polar cap and not to the absolute value of the 
scalar potential, as had been generally assumed. This provides a 
wide range of opportunities for comparing the observational data 
on pulsar profiles with the results presented here and in our pre- 
vious investigations of oscillating magnetopsheres. Stellar oscil- 
lations, indeed, affect the scalar potential in the polar cap region 
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and are expected to have a significant effect on the character of 
the subpulse drift. We plan to devote our future research to a 
closer analysis of the observational aspects of these ideas. 

Acknowledgements. V.Morozova would like to acknowledge the support of the 
German Academic Exchange Service DAAD for supporting her stay at ZARM. 

Appendix A: Extended geometrized system of units 

We recall that the definition of geometric units of time and 
lengths is obtained by setting the speed of light c and the grav- 
itational constant G to pure numbers. This implies that seconds 
and grams of the cgs system can be written as 



Is = 2.997924x10' 



- cm 

c 



Ig = 7.424157x10"-'' ( — 



(A.l) 



(A.2) 



Within this general setup, a convenient unit of space is still re- 
quired. The cm is of course a bad choice in a gravitational frame- 
work and the gravitational radius = GM/c^ is instead chosen. 
In order for this new unit to be more convenient than the cen- 
timeter, the mass M of the system has to be sufficiently high. 
From the physical value of the solar mass and Eq. ( IA.2b . we find 
the relations between the cgs units and the new unit of length rg 



icm = 6.772289x10"*' rg , 



Is = 2.030281x10" 



Ig = 5.027854x10 



-34 



(A.3) 
(A.4) 

(A.5) 



In the traditional geometrized system c and G are set equal to 
unity. The equations Eqs. (IA.3b - (IA.5l l allow us to perform the ge- 
ometrization of any dynamical quantities. Quantities involving 
the charge, on the other hand, are naturally geometrized within 
the Gauss system of units, which we adopt here in addition to 
the choice c - G = 1 . In Gauss units, the unit of measure of the 
magnetic induction, the Gauss G, is given by 



-1/2^-1 



IG = g'/^cm 
which, by virtue of ( IA.3b -( lA3] ), provides 



fi„eo =4.2439x 10--" -^IBg, 



M 

Mo 



(A.6) 



(A.7) 



where Be and Bgeo are the pure numbers expressing the intensity 
of the magnetic field in Gauss units and the geometrized system, 
respectively. We note that in these units the charge-to-mass ratio 
is a pure number, and, in the case of electrons, is given by 
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